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ABSTRACT 


Experimental results are presented for saturated 
nucleate pool boiling of distilled water and ethyl-alcohol 
from 1-7/8" diameter Nickel disks kh pool depths from 1.0 
inch to a thin film of 0.015 inch. Experimental runs were 
conducted on a mirrorlike surface utilizing distilled water 
and ethyl-alcohol; on a grooved surface with distilled water 
as the working fluid; and on a surface covered wiheaeo  lavews 
of Nickel 50 x 40 mesh with distilled water as the working 
EEC . a 

Results show that liquid level has little effect on the 
heat transfer coefficient above a level of 0.2 inch on all 
three test surfaces; however, below 0.2 inch a rapid increase 
in the heat transfer coefficient occurs with the lowering 
of the liquid level. Two types of thin film dryout were 
noted on the mirrorlike surface: dryout under a dome and 
dryout in absence of a dome. A grooved surface appears to 
provide a means of obtaining a stable high heat transfer 


coefficient at liquid films as thin as 0.015 inch. 
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lL. INPRODUCT ES 


A. BACKGROUND 

Since the advent of the space age numerous investigations 
into the theory and applications of heat pipes have been con~ 
ducted. Although noted by Gaugler in 1942 [1], the relatively 
simple principle of heat pipes did not come into its own 
right until 1964 with the studies of Grover, et al. eae 
In general, a heat pipe consists of a closed eva. ima Tr lead: 
container in which vaporization and condensation of a fluid 
takes place. A porous wick material lines the container 
while maintaining a vapor space in the center. Heat added 
at one end of the container causes evaporation of saturated 
liquid in the wick. The vapors then travel through the 
center space to the heat sink where the vapor condenses aa 
the wick. The resulting condensate is returned to the heated 
end of the container by the action of capillary forces in the 
liguid layer and aie cycle continues. Heat pipes are poten- 
tially inexpensive heat transfer devices which can operate 
at high heat fluxes, have no moving parts, are independent 
Guegravity, require no nee, ona operate silently and 
reliably for extended periods of time. 

Garg [3] conducted an investigation of heat pipe technology 
‘for Naval applications. He concludes that heat pipes may be 
useful for a wide range of Naval applications because oteenee 
properties cited above. Also contributing to Naval application 


are their compactness and light weight for a given use when 


a 





compared to conventional heat dissipation devices; their 
flexibility in configuration, which can be used to transfer 
heat to or from otherwise inaccessible locations; and their 
ability to maintain a constant temperature over a large 
Surface at varying power levels. At the present time, space 
application remains the main area in which heat pipes are 
being utilized; however, Garg points out that several of the 
Naval laboratories are looking at the heat pipe for utility 
yee applications although no work in thiswarea *has»commenced. 
A couple of the uses suggested by Garg for Naval application 
are: 


1) Heat pipes could be used to protect the 
refractory walls of boilers. 


2) Cooling of large power transformers could 
be carried out more effectively by using heat 
Pipes in place of the presently used radiators. | 
3) Heat pipes could be efficiently used where heat 
is to be removed from or added to high voltage 
or high radiation areas or spaces. 

Heat pipes have been operated using various temperatures, 
pressures, wicking materials and fluids [4]. Studies have 
Shown that the heat pipe performance is mainly limited by the 
wick resistance of the capillary structure when uSing ordinary 
fluids. When the capillary action is exceeded the supply of 
working fluid to the evaporator is insufficient and dryout 
occurs in the evaporator causing failure of the,pipe. Nucleate 
boiling may also occur in the wick structure leading perhaps, 


fO GEYOuL. 
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The mechanisms of evaporator heat transfer in a heat 
pipe have received considerable attention. Alleavitch [5] 
designed an apparatus which would simulate the vaporization 
zone of a heat pipe in order to study this heat transfer 
mechanism. Ferrell and Johnson [6] reported their results of 
an experimental investigation of the mechanism of heat trans- 
fer during the evaporation of a fluid in a porous wick 
structure in contact with a heated surface. Seban and 
Abhat [7] investigated evaporator heat transfer from screen 
wicks. Marto and Mosteller [8] studied the effect of nucleate 
boiling on the operation of low temperature heat pipes, and 
its effect on the mechanism of evaporator heat transfer. 
Gregory [9] investigated nucleate boiling from a mesh covered 
surface simulating a heat pipe. 

Since heat pipes operate at low liquid levels, it is 
appropriate that an investigation of nucleate boiling at low 
liquid levels be made. Previous attention toward this area 
of nucleate boiling has been small, with the majority of 
studies on nucleate boiling being made from a surface situated 
well below the liquid surface. Nishikawa, et al. [10] in- 
vestigated nucleate boiling at liguid levels from 1 mm to 
meemm on Several flat horizontal surfaces using various 
working fluids. They suggest that a critical layer thickness 
exists below which dryout will always occur. Patten and 
Turmeau [11] conducted heat transfer experiments from a flat 


horizontal surface covered by liquid layers from 5 mm to 0.5 mm. 


lig 





and proposed mechanisms for nucleate boiling in thin 


liquid layers which differ from normal pool boiling. 


Be LHESIS OBJECTIVE 

The objective of this study was to systematically 
investigate nucleate boiling at low liquid levels ona flat 
horizontal surface, and to determine what effect test surface 


Penditions have upon thin film boiling. 
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Ii<«. EXPERIMENTAL S DES EG 


fiw FACTORS CONSIDERED AND CHOICE OF MATERIALS 

The overall design of the boiling apparatus used in this 
study was largely influenced by the requirements to allow 
the condensate to return directly into the boiling area, and 
to allow close monitoring of the liquid level in the boiling 
area. While striving to satisfy the above requirements, 
other factors were kept constantly in mind. Some of these 
factors included: permitting visual observation of the boil- 
ing process, flexibility of operation, simplicity in design 
and Bee ructicn, as well as safe operation of the system. 

The final design was quite similar to that used by 
Gregory [9]. In fact, several of his components or design 
ideas were used with some modifications. The final design 
consisted of two major components: a test surface/heater 
assembly which contained a nickel cylinder to boil from and 
the necessary heater elements to provide the’ required heat 
fluxes; a glass condenser assembly which provided a medium 
Peenold the test fluid and to permit the passage of external 
coolants. These components will be described in detail in 
the following section. 

Previous studies of heat pipe operation at the Naval 
meotgraduatesSchool [9] and several StudWes dwring the period 
of this report, had used nickel pipe with nickel screen 
serving as the wicking material and for this reason, a 


nickel boiler surface was used. Nickel screen, fed by gravity, 
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was used as the wicking material in some runs because of its 
ease of fabrication, dimensional stability tremmaiecongue. 
tivity, and wetting properties as recommended by Kunz, et 
cee {12}. 

With the exception of the nickel and electric resistance 
heaters, all materials used in this experiment were readily 
available and were used for simplicity of construction as 
well as to avoid risk of corrosion and contamination of the 
boiling fluids - distilled water and ethyl alcohol. Figure 1 
is a cross-sectional view of the boiling apparatus, and 


Figure 2 an overall view of the experimental apparatus. 


B. COMPONENTS 

The two nickel cylinders used by Gregory [9] were 
modified for use in this experiment. Originally the nickel 
cylinder was machined from a 2 1/2 inch Nickel 200 rod toa 
iors inch diameter with a 0.1875 inch Vip, 0.050 inch thick, 
extending from the sides of the cylinder, 0.010 inch from 
the boiling surface itself. The lip was undercut until the 
final thickness at the machined diameter of the cylinder was 
eeo'> inches. “Modification included the leveling of the boil- 
ing surface to the lip and the shortening of the overall 
ferwgit to 2" from 2 1/4" with the bottom 0.25 inch threaded: 
The bottom of the lip provided a means of supporting the 
heater-boiling surface assembly from the phenolic base while 
the undercutting reduced the extent of conduction losses 


through the lip. 
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Six thermocouple wells, 0.041 inch diameter, were 
Grilled at a number of axial and radial positions in the 
nickel cylinder to accommodate the thermocouples. Four 
of these wells were drilled at ninety degree angles to each 
other, .500 inch below the boiling surface to a depth of 
-625". The fifth hole was drilled 1.000 inch below the 
surface and the last hole was drilled 1.500 inch below the 
surface, both to a depth of .625". No two holes lay in the 
same plane parallel with the longitudinal axis of the cylinder. 
The first four were placed a plane perpendicular to the 
axis to determine if the temperature distribution throughout 
the nickel was uniform. 

As a support for the nickel cylinder and as a base for 
the boiling apparatus a 8.5 inch phenolic disk 5/8" thick was 
chosen for its low thermal conductivity, low porosity, and 
lei machinability. The 1.875 inch hole in the center was 
Seumeerbored 0.050 inch to a radius of 1.125 inch so that 
the actual boiling surface would sit flush with the phenolic 
surface and was counterbored on the bottom .250 inch to a 
radius of 1.562" to accommodate the thermocouples. A place 
for the nickel screen ring was counterbored .125 inch from a 
mezo0' ~- 1.437" radius*om thesphenolic dish. Groeéwes* weue 
cut .1875" in the disk to accommodate both the inner and outer 
glass fluid containers along with their '0O' ring sealant. 

The 3-inch high copper heater unit, chosen for its high 


thermal conductivity, was machined to a diameter of 1.875 
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maehes, threaded at the top and had seven 023575 inehehove>s 
2.250 inches deep, drilled into the bottom. Watlow Firerod 
Cartridge Heaters, 0.375 inch in diameter, 2.250 inches long 
and having a capacity of 250 watts each, were inserted into 
the holes. The bottom surface of the nickel cylinder and the 
top surfaces of the copper were polished until uniform, then 
coated with Silver Goop to ensure good contact and to prevent 
sticking after heating. A threaded nickel ring held the two 
cylinders in direct contact. A stainless steel collar was 
attached to this nickel ring to which 6 hold down bolts were 
attached to form a sealant between the nickel boiling surface 
and the phenolic base. (Figure 1) 

After assembly, the cylinders were covered with 1l-inch 
thick asbestos preformed insulation which was covered by 
several layers of Flex - Min K asbestos insulation. The 
Gemwer portion of the apparatus, i.e. cylinders and insulation, 
was contained ina 61/4" diameter stainless steel container 
meunmewed on top of a leveling device which consisted of a 
bottom plate with 3 jacking screws attached, a l-inch high 
stand and a hold down plate for bolting down the nickel 
cylinders and copper heating assembly. Two level indicators 
were mounted on the phenolic support for the nickel cylinder 
which, along with the jacking screws, provided a means of 
leveling the boiling surface. 

The condenser assembly consisted of inner and outer fluid 


containers cut from 6" high glass tubes; 4.140 inches in 
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diameter by .110 inches thick, and 5.290 inches in diameter 
by .130 inches thick respectively. An aluminum top cap, 

mm cin diameter, 7/16" thickswas, placed on @topwor tie 
glass condenser assembly. A rubber gasket formed a sealant 
at the top and bottom of the assembly. If a coolant was used, 
the fluid entered the glass condenser through a hole in the 
phenolic base and left through a hole in the aluminum cap. 

The boiling fluid entered the boiling area through a hole in 
the aluminum cap and could be removed through this hole with 
an eyedropper during the constant heat flux experimental runs. 
Six bolts provided a means of sealing the top cap, glass con- 


denser assembly, and phenolic disk together. 


C. INSTRUMENTATION 

Temperature measurements were made with copper-constantan 
thermocouples using a distilled water-ice bath as a reference. 
A total of six thermocouples were used to measure the tempera- 
ture profile in the nickel cylinder, while one thermocouple 
was used to monitor the fluid temperature. Thermocouple out-, 
put was measured with a Hewlett-Packard 2010C Data Acquisition 
, System containing an integrating digital voltmeter and guarded 
data amplifier. Thermocouple selection was provided within 
the DA system. This instrument gave a printed digital output 
with an accuracy of + 0.5 MiGhOvo res. 

Heat flux control was provided initially by means of a 


Wariaac to regulate the voltage input alenguwich a vole. eer 
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and ammeter to monitor power input. A D.C. regulated power 
supply was later used to provide this control. 

A micrometer fluid level indicator was mounted on the 
aluminum cap which permitted measuring of the liguid level 
feeectly over the center of the boiling surface and radially 
outward to 1.125 inches. This device consisted of two micro- 
meters mounted in an aluminum block. One micrometer, mounted 
horizontally, measured the radial distance while the other 
micrometer, mounted vertically, measured the liguid level 
with its probe. An electrical circuit between the vertical 
micrometer and test surface would close and activate a light 
when the probe made contact with the test fluid, signalling 
the point to read the micrometer which measured the liquid 


level. 
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Jil. EXPERIMEN@AL PROCEDURE 


Poo PREPARATION OF TEST .SURBPACES 

Throughout this research work the condition of the test 
surface had to be Known in order to compare results for re- 
producibility. No theoretical nor experimental method exists 
today which can predict the number of active sites of a given 
size on a test surface which has undergone a prescribed 
polishing procedure. However it is hoped that the structure 
of a given test surface may be adequately defined if an 
accurate description is given of the polishing procedure used 
in preparing the surface [13]. Much time was consumed in 
preparing the boiling surface as explained in the following 
paragraphs. 

tee Mirrorlike Finish 

In the majority of tests a mirrorlike finish was 
required. In order to obtain this finish, great care had to 
be maintained at all times. Rough polishing was done on a 
Buehler Metallurgical Manual Grinder. This grinder contained 
strips of 0 emery polishing paper which were placed over 
rectangular glass plates and held in place by clamps. The 
test surface was placed face down on the 0 emery paper and 
was stroked over the paper in one direction Only, White) vada 
the original scratches were removed. This required approx- 
—t 

imately 200 strokes and several changes of the polishing 
paper. The test surface was then rotated 90 degrees and was 


carefully stroked in this direction over the 2/0 emery 
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polishing paper until all the scratehes left Emom Ehen0) emenm, 
polishing paper were removed. Upon completion of the rough 
polishing, the surface had a high gloss finish although there 
were still many visible scratches. 

The test surface was then cleaned in an ultrasonic 
cleaner, which helps to wash out small dust particles and 
bits of metal trapped in the cavities on the surface. This 
was followed by washing with an alcohol jet and drying with 
a hot air jet. It was then put on the No. 6 diamond compound 
wheel (6 microns). It is important that the diamond compound 
be put on the center area of the grinding wheel and be diluted 
myealcohol before starting the grinding operation. The test 
surface was rotated 90 degrees and was held gently near the 
edge area of the wheel. Alcohol was constantly injected on 
the wheel cloth. Operation was continued until the scratches 
due to the 2/0 emery polishing paper disappeared completely. 
The test surface was then cleaned again, and put on the No. 3 
diamond compound wheel (3 microns) with a 90 degree rotation. 
The process continued to a No. 1 diamond compound (1 micron) 
wheel. During the last few minutes of grinding process on 
the No. 1 wheel, the test surface was placed near the center 
part of the wheel where the rubbing speed is lower, and a 
heavier pressure was applied. After a few minutes the surface 
would become shining, mirrorlike smooth. It was once more 
Washed in an ultrasonic cleaner, alcohol jet, and@heteain jee. 
The surface was then carefully covered with masking tape until 
meady fOr use. 
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If refinishing was desired after an experimental 
run(s) the above process was repeated starting with the #3 
diamond wheel, and continued through the #1 diamond wheel. 
Table 1 indicates when refinishing was accomplished prior 
ieona run. 

2. Grooved Surface 

EXperasie Wat eUnic wee Ou OG NGccmamsed grooves in 
the test surface. This was accomplished by machining the 
same nickel cylinder as used in runs iMthGewaee/ with a 
.008" mill saw, to a depth of .006" with a nominal 62 grooves 
per inch. Figure 3 shows a top view and side view of the 
finished test surface. After machining, the surface was 
cleaned in an ultrasonic cleaner, followed by rinsing in 
distilled water, alcohol, and reagent grade acetone and 
finally dried in a hot-air jet. 

3. Mesh Covered Surface 

Experimental runs 36 and 37 required wire mesh to be 
Sintered to the pase Surface. Four layers of nickel, 40 x 50 
wire mesh screen, were Sintered to the surface of the nickel 
test surface used in runs 28 - 31, in a Watkins Johnson 
hydrogen controlled atmosphere furnace for two hours at 
1200°C. The wire size was 0.009 inches, with a percent open 
area of 35.3. Thickness of the 4 layers after sintering was 
.081 inches. 

Prior to sintering, the.wire mesh was prepared as 


follows. The pieces were cut to a diameter of 2 1/4 inches 
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and cleaned in an ultrasonic cleaner, followed by rinsing in 
distilled water, alcohol, and reagent grade acetone. This 
prepared the samples for thermal oxidation as recommended by 
Kunz, et al. [12]. Thermal oxidation was accomplished by 
placing the screen pieces in a Lindberg Hevi-duty Electric 
Furnace for two hours at a temperature of 500°F. Thermal 
oxidation results in the formation of nickel oxide on the 
screen which, having a smaller contact angle than pure 
nickel, increases the wettability of the nickel screen. 

Since the lip of the test surface used in the sinter- 
ing process had previously been damaged, it was decided to 
Sinter the mesh over the center 1.8 inch diameter of the test 
surface, and allow the rest of the 2 1/4 diameter screen 
pieces to rest on the lip area only. A special stenetiy 
device was made from an asbestos ceramic material and Nichrome 
Wire, which applied pressure over the sintering surface. Care 
was taken to orient the screen relative to the mesh opening; 
however, during one Sintering process the mesh slipped, re- 
sulting in no specific screen orientation during the runs 
with the sintered surface. Other than washing the surface 
with acetone, no special preparation was done on the nickel 


surface prior to sintering the screen to it. 


B. PREPARATION OF EQUIPMENT 
Once the test surface was prepared, assembly of the 
boiling apparatus proceeded. All parts which came into con- 


tact with the test surface or test fluid were thoroughly 
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cleaned and degreased. This process consisted of using 
Trichloroethylene for degreasing, followed by a detergent 
washing, a distilled water rinse and a final rinse with 
reagent grade acetone. 

The assembly process consisted of assembling the lower 
section of the boiling apparatus followed by a check of the 
voltmeter and ammeter for proper connections, operation and 
zero readings. The thermocouples in the nickel cylinder 
were then checked by applying an input of 25 volts and read- 
ing the Data Acquisition System. Throughout the data 
collection process, the Data Acquisition System was left 
operating to prevent any malfunction caused by repeated start- 
up and securing. As a result, only a minor calibration, 
zeroing of the apparatus, was required prior to each run. If 
the instrumentation checked out, the top section of the boil- 
ing apparatus was assembled and checked for leaks. Any 
external leaks were corrected using a silicone rubber sealant, 
while internal leaks required disassembly of the top section 
to investigate the cause. Changing of the test surface 
required the above procedure to be repeated. Prior to the 
beginning of any experimental run the leveling device was 


@iiecked to insure a level test surface. 


C. NORMAL OPERATIONS 


ih =Constamt Liquid Level 


The system was prepared for experimental runs by 


gntroducing sufficient test fluid into the testearcasee 
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bring the level up to the desired value. The level in the 
test area was monitored by the liquid level device mounted 
‘on the top aluminum cap for levels below 2" and by means of 
a ruler for levels over 2". A heat flux corresponding to a 
40 volt input (approximately 180 watts) was added and the 
test fluid was allowed to boil for 30 minutes in order to 
degas the fluid. Power was then secured and the system 
allowed to cool until the thermocouple just above the heaters 
read less than 3.200 millivolts (approximately 170°F). A 
final check of the fluid level was made and the system was 
considered ready for operation. Runs 1-3 used tap water as 
a coolant, while during run 4 an air blower was used to pro- 
vide coolant for condensing the vapor. 

At the start of each run, the atmospheric pressure 
was recorded from an aneroid barometer and used to calculate 
the saturation temperature of the test fluid. Voltage was 
increased in 5-volt increments starting with 25 volts. Data 
was collected at the 5-volt levels, up to a maximum of 80 
volts (approximately 500 watts). At this voltage, the fluid 
level fluctuated so rapidly, that the system couldn't condense 
most of the vapor. Data consisted of 6 thermocouple readings 
in the nickel cylinder, 1 thermocouple reading in the liquid, 
l thermocouple reading in the cooling chamber, the input 
voltage and input current. Data was obtained when a semi- 
steady state SeeEat ion for a given heat flux was reached. 


Semi-steady state was assumed to be the point at which the 
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wOoltage of the thenmocouplesneamesteene Vicarte go awaoee 
creaSing at a rate less than .010 millivolts per minute, which 
corresponds to a temperature rise rate of approximately .3°F 
per minute. Between each data point the water level was 
measured. 

Throughout the experimentation the apparatus was 
Monitored closely and significant occurrences and visual 
observations were noted on the data sheet. 

yee COnNStant Heat Flux ; 

Problems with maintaining the fluid level were en- 
Gountered at high heat flux during constant liquid level runs 
requiring the experimental procedure to be modified. Instead 
of varying the heat flux while maintaining a fixed fluid 
level, the heat flux was held constant by using a 40 volt DC 
regulated power supply, and the test fluid level was varied. 
Fluid levels from 1.00" down to .015" were used in these 
tests, while constant heat fluxes from 8,000 - 25,000 
Btu/Hr-Ft¢ were applied. Power to the heater block was ad- 
jJusted to give the desired heat flux. Steady state conditions 
were assumed to have been met when there were virtually no 
meeetuations (+ .003 MV) in the thermnocouple millivolt 
readings for a period of 30 minutes . Thermocouple readings 
and the liquid level readings were then recorded. Because 
condensation formed on the liquid level probe and vapor domes 
were present on the liquid surface, three readings of the 


liquid level were taken and the average value recorded on 
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the data sheet. The liquid level was then reduced to the next 
desired value, and after conditions stabilized, thermocouple 
readings and liguid level measurements were again recorded. 
This process was repeated until dryout occurred on the 

boiling surface. The power to the heater block was then 
adjusted to give a new desired heat flux, and more test fluid 
was added to bring the liquid level to approximately 1.00" 


and the above process repeated for this heat flux. ms 


D. DATA REDUCTION 

The thermocouple millivolt readings, applied current and 
voltage readings, and atmospheric pressure readings were put 
into a data reduction program written for an IBM 360 computer. 
This program calculated the thermocouple temperatures, the 
heat flux from the power input, and thermocouple temperatures, 
Mme conductivity of the nickel cylinder, the BEE IOn 
temperature, the surface temperature, surface-solution 
temperature difference as well as plotting graphs of heat 
flux vs. surface-saturation temperature. Because of the 
linearity of the temperature distribution in the nickel 
@ytinder the heat flux could be calculated by using Fourier's 
heat conduction equation: 
Swan (1) 

A L\X 

The surface temperature was determined by extrapolating this 
itineal temperature distribution to the test suriace as PS nileubncel 


wn Appendix D. 


Zo 





A sample calculation of the data iS given in Appendix B, 


along with an uncertainty analysis in Appendix C. 
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Iv. RESULTS AND DISCUSSION 


In all, a total of 37 experimental runs were completed. 
Table 1 presents a summary of all experimental runs. Results 
have been divided into 6 categories for discussion purposes: 

1) Testing of Experimental Apparatus 
2). Kspect Of Boiling 
3) Reproducibility and Aging Effect 
4) Effect of Liquid Level on Heat Transfer Coefficient 
5) Vaporization of Liquid Levels Below .200 Inch 
6) Special Surfaces to Improve the Heat Transfer 
Coefficient 
A. TESTING OF EXPERIMENTAL APPARATUS 

The purpose of Run 1 was to test out the design apparatus 
emda to check the data reduction program. It initially ap- 
peared that the apparatus was functioning correctly and that 
the results obtained were similar to previous studies; however, 
by the end of Run 3 it was Seedient that the use of the tap 
water coolant for condensing purposes prevented constant 
Peol boiling and interfered with the normal pool boiling 
curve. Figure 4 shows the results of Run 1 and 3. At high 
heat flux during Run 3 the boiling water circulated ina 
system, flowing upward in the central part and then flowing 
downward along the container walls transferring heat in the 
process. This sub-cooled water flowed inwards toward the 
test surface which resulted in the change in the normal 


boiling curve. This problem could have been corrected 
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Heat Flux X 1073, Btu/Hr - Ft 
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utilizing an enclosed heat exchanger-pumping system capable 
of maintaining the coolant water at a few Blegon Se less than 
the saturation temperature; however, this equipment was not 
readily available and an alternate method was tried. Results 
of Run 4 (Figure 4) showed an improved boiling curve when 
forced air was used as a coolant; nevertheless, a decision 
was made to eliminate all forced coolant and leave the glass 
exposed to the ambient temperature which appeared to be 
sufficient to condense the vapors. 

The heat flux computed from the power input was compared 
to the heat flux based on the temperature profile in the 
nickel cylinder for Runs 4 and 5. The results indicated a 
heat loss of 13% and 19% for these runs respectively. The 
Majority of the heat loss was through the insulation and thus 
_— a function of the manner used in packing the insulating 
material which varied every time the apparatus was taken apart. 
Because of this substantial power loss, the heat flux calcu- 
Jated from the re eeee gradient was the one used in 
correlating all subsequent data. 

The heat transfer apparatus and experimental techniques 
were then evaluated by obtaining data on pool boiling of 
distilled water. Figure 5 shows pool boiling data for Run 5 
at .995" compared to Rohsenow's correlation for pool boiling 
with a surface-fluid combination of water-nickel (C.;, = .006) 
[14]; to the pool boiling data of Kunz, et al. [12] from an 
inconel flat plate; to the pool boiling data from a square 


Stainless steel surface of Alleavitch[{5]; and to the pool 
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Bowling data of Gregory [9] utilizing a similagececcmaeoac noe 
The displacement of the data from Rohsenow's correlation was 
assumed to be due to the difference in the nature of the 
working fluid-surface interface. The displacement of the 

data from those of Gregory is believed to be caused by the 
difference in surface preparation. Gregory's surface was 
roughened with No. 240 Three-M-Ite Elek-Tro-Art cloth to 
promote nucleate boiling, while the metal surface for this 
study had been polished on a diamond wheel to 1 micron. 

Figure 6 shows the effect of pool depth on the boiling 
G@meve. Pool depths of .995, .495, and .255 inches were 
@empared. The curve shift to the right at .255 inches indi- 
cates a substantial decrease in the heat transfer coefficient 
at this liquid level; however, the cause of this shift was 
not known. Rohsenow [14] suggests that an aging effect may 
exist even though care was taken in repolishing the surface 
using the same repolishing technique as in Runs 5 and 6. 

| As the heat Pee increased it was impossible to keep the 
water level constant, so another method was €ried ian whiched 
constant heat eles was applied and the water level varied. 
This was similar to the methods used by Patten and Turmeau [11] 
and Nishikawa et al. {10]. Runs 8 and 9 were used to evaluate 
this method at two heat fluxes. It was noted that the A. C. 
input power fluctuated randomly anaes D.C. regulated power 
Supply was acquired in order to stabilize the power input 


fOr futwre runs. 
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Bee ASPECT OF BOILING 

During normal pool boiling, water circulates ina 
specific pattern; flowing upward in the center area, downward 
along the container walls and finally flowing inward toward 
the center of the heating surface to repeat the process. This 
type of pattern was found to occur in levels above 1 inch as 
shown in Figure 7 for a water level of 2 inches, and is 
characterized by the pyramid form the flow of bubbles project. 
Below 1 inch this specific pattern ceases and the bubbles 
climb straight up to the surface, as shown in Figure 8. 

mm the liguad range, from 1.0% to 02:2 inehes, sit was nored 
that vapor domes formed on the liquid. These domes formed 
as the result of the coalescing of small bubbles as they 
leave the boiling surface. Vapor dome formation appeared to 
be influenced by two factors, the liquid level and the heat 
flux. At a constant heat flux the Size of the dome increased 
as the liquid level was decreased. This increase in dome 
size appears to be a result of the decrease in liquid circu- 
lation as the level was decreased. While at a constant liquid 
level, dome size increased with increasing heat fluxes as a 
result of the increased number of nucleating sites. Figures 
9 and 10 are photographs of a vapor dome forming at a liquid 
evel of 1/2 inch and a heat flux of approximately 15,000 
Btu/Hr - Ft2. These domes are l 1/4 - 11/2 inches in diameter. 
In several cases, one vapor dome would cover the entire heating 


Surface. These extremely large domes were noted at heat 
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fluxes of 18,000 - 25,000 Btu/Hr - Ft? and at liquid levels 
between .2 ~- .4 inches. Stability of these domes decreased 
as heat flux increased; however, in the .6 - .2 inch liquid 
Beevel with a constant heat flux, the stability remained 
constant with decrease of liquid level. Domes remained in 
foeteror 10 - 15 seconds at lower heat fluxes and l - 2 
seconds at higher heat fluxes used in this study. 

Since bubbles form from nucleating sites, an attempt was 
made to investigate the effect of heat flux and liquid level 
on the number of nucleating sites. It was quite evident from 
ene runs at a constant liquid level that as heat flux is 
increased, the number of active nucleating sites increased. 
This increase in nucleating sites was not evident however 
under conditions of consizant headke Tt lux while lowering the 
liquid level. 

When ethyl-alcohol was used on the test fluid, ee 
nucleating sites appeared on the rim of the test surface and 
formed a 360° curtain around the active sites on the surface, 
thus interfering with any flow pattern the bubbles might have 
hemmed. It Pes noted that the bubbles released from rogue 
emees Climbed straight up to the surface at liquid levels 
below l inch. Levels above 1 inch were not checked. 

With ethyl-alcohol no large vapor domes formed. In fact, 
domes larger than 1/4 inch were rare. The population of domes 
mas much greater Chan with distilled water beeduse of Gem 


coalescing of bubbles in water. Individual bubble diameters 
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ee 


appeared to be approximately 1/3 the size of bubbles formed 
in distilled water. The increase in nucleating sites and 
bubble population as a result of an increase in heat flux 


[ere much more evident than in distilled water. 


©. REPRODUCIBILITY AND AGING 

The time required to disassemble the boiling apparatus, 
remove the test surface and refinish it in order to reproduce 
results, waS approximately 8 hours. This time, although well 
spent, was considered excessive and another method for repro- 
ducing results was desired. Rohsenow [14] points out that 
experience teaches that aged surfaces have higher required 
AT for a given heat flux. Initial runs carried out in this 
study under Similar conditions indicated the presence of 
Meeemaging effect on the test surface as shown in Figure I]. 
It should be noted pa during these aging runs no attempt was 
made to clean or refinish the test surface and that distilled 
water was used as ing test fluid. This aging effect may be 
caused by a scale or deposit formation from the boiling 
liquid, or as a result of oxidation or other chemical reactions. 
When this happens the nucleating site(s) may shrink requiring 
higher superheat for activation. 

Additional runs were conducted to determine if this 
aging effect could be minimized after several hours of boiling. 
Results shown in Figures 12 and 13 indicate that the aging 


effect diminishes with boiling time. These results also 
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indicate that when the aging effect subsides at one heat flux, 
it appears to subside at all heat fluxes. Higher heat fluxes 
appear to age the surface much faster, however. 

Thus reproducibility of results was found to be possible 
once the surface had aged. Reproducibility between two 
surfaces was attempted. Two nickel surfaces were machined 
similarly and polished utilizing the same technique. Results 
of tests on surface 1 were discussed above, BSL BITTE 2 was 
tested to verify these results. In assembling the boiling 
apparatus the rim of surface No. 2 was acclagenecly bent orcs 
sulting in a dish-shaped surface profile. The appearance one 
cracks at the rim/surface interface may have acted as 
additional nucleate sites which were not present on surface 
#1>- therefore, no attempt was made to compare results to 
see if reproducibility between surfaces was possible. How- 
ever, 4 runs were conducted on this surface to verify the 
aging effect noted on surface “41 and the ability to reproduce 
results on the same surface. As shown in Figure 14, the 
same aging effect and ability to reproduce results on the 


same surface occurred on surface 2. 


nf 


me EFFECT OF LIQUID LEVEL ON HEAT TRANSFER COEFFICIENT 

An investigation into the effect of liguid level on the 
Heat transfer coefficient was conducted along with the study 
of reproducibility and aging. Results are shown in Figure 15, 
which is a graph of the heat transfer coefficient (h, Btu/Hr - 
Ft? - °F) vs. liquid level (inches) élbbaabcley IOUS ee OS constant 


heat fluxes on a surface where the aging effect has subsided. 
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In actuality, the constant heat flux was assumed to be the 
everage heat flux, determined from Ehe Genpemeeurese.erile, 
encountered during a run, which resulted in an uncertainty 
Poet 400 Btu/Hr = Ft? Liquid levels betweem 1 inch and 
0.040 of an inch were used in this study. 

bigquid levels in the range 0)..2 - 1.0 anem appear tovhave 
little effect on the heat transfer coefficient at the heat 
fluxes studied. In this level range the heat transfer co- 
efficient would vary about the same average walue, approximately 
Eeoebcu/Hr — Ft? °F. Of course, aS expected, the higher the 
heat flux, the greater the heat transfer coefficient. One of 
the interesting points noted in Figure 15, amd also noted in 
Figures ll, 12, and 13, was the fact that the heat transfer 
Seetficient decreased slightly in the range 1.0 - 0.75 inches; 
then increased slightly to a maximum in the 0.4 - 0.75 inch 
range; followed by a decrease to a new low heat transfer co- 
efficient in the 0.2 - 0.4 inch range. As noted in Run 3l, 
as shown in Figure 14, the above heat transfer seam Cue 
results occurred on a second test surface under an identical 
Situation. This appears to be in partial agreement with the 
results of Nishikawa, et al. [10] except in the 0.4 - 0.75 
inch liguid level, where they had the heat Matenstier coeffi- 
cient decreasing by slow degrees with decreasing liquid level 
an the 1.2 inch - 0.4 inch range. No firm conclusions can be 
drawn here until additional research is completed. 

Levels below 0.2 appear to have a tremendous effect upon 
the heat transfer coefficient as shown in Boye@re 15. An 


abrupt increase in the coefficient was noted, in most cases, 


>] 





below the 0.2 inch liquid level. In each case where this 
increase was noted, dryout occurred on the test surface. 
Dryout was characterized by a sudden fluctuation in the 
thermocouple temperatures, and by a bright spot appearing on 
the surface. In this range of liquid levels the lower the 
heat flux the more remarkably the liquid level affects the 
heat transfer coefficient. It should be noted here, that 

in Figure 15 it appears that dryout initially occurred at 
the same liquid level for all heat fluxes. This is strictly 
a cOincidence because no attempt was made to find the exact 
level where dryout commences and the uncertainty in the 
iieuiad Level at a heat flux of 20,000 Btu/Hr - Ft? was 

fee o> inch. 

Figure 16 shows the effect of pool boiling at liquid 
iewels Of 1 inch, 1/2 inch and 1/8 inch. The pool boiling 
curve at 1/8 inch shows an abrupt shift to the left indi- 
Cating a much higher heat transfer coefficient. This coin- 
cides with the ee of Gregory [9] and Nishikawa, et al. 
[10]. | 

The effect of liquid level on the heat transfer coefficient 
was also investigated with ethyl-alcohol as the test fluid. 
Five runs were conducted at various constant heat fluxes 
with the results shown in Figure 17. The results were similar 
to those described above with distilled water but with much 
lower values of heat transfer coefficient, as expected. The 
effect of liquid levels in the 1.0 - 0.4 inch range appears 


te be rather small. No increase in Heat trancter coefficient 
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was noted in this 0.4 - 1.0 anch liquid levels tise reecrcun. 


with the results of Nishikawa, et al. [10]. A slight decrease 
occurred in the 0.4 - 0.2 inch range with an abrupt increase 
mathe coefficzent occurring in the 0.080 - 0.120 inch range. 


This increase in heat transfer coefficient was accompanied by 
dryout of the test surface in most cases; however, this was 
much harder to see than in the distilled water case because 
no large vapor domes formed, and the large number of nucleat- 
ing sites present prevented good visualization of the boiling 


surface. 


E. VAPORIZATION OF LIQUID LEVELS BELOW 0.2 INCH 

As mentioned before, liguid levels below 0.2 inch appear 
to have a tremendous effect upon the heat transfer coefficient. 
With water as the test fluid, dryout coincided with an abrupt 
increase in the heat transfer coefficient. Dryout occurred 
in two ways: dryout under a dome and dryout in the absence 
of a dome. The first way, dryout under a dome, is the result 
of internal pressure developing within the dome forcing fluid 
meow the dome aside until the test surface 4S exposed under 
Eae dome. tn addition the presence of the dome prevented 
bubble microconvection effects across the test surface. 
This type of dryout was present throughout the liquid levels 
Studied below .200 inch. It was characterized by a sudden 
fluctuation in the thermocouple temperatures, and a bright 
spot appearing under the dome. At levels above 0.060 inch, 
even with this type of dryout occurring, it was possible to 


Sustain nucleate boiling. Figure 18 shows vapor domes 
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forming at 0.090 inch of water, with dryout occurring under 
some domes. 

Below 0.060 inch, the second type of dryout, Grayour 
without a dome, would SOCGE along with dryout under a vapor 
dome. This type of dryout was characterized by a much higher 
drop in surface temperature along with artastes spreading 
rate. This spreading could be controlled for a while in the 
range of 0.040 - 0.060 inches by means of the vapor domes 
pushing fluid over the exposed surface; however, around 0.040 
inch dryout without a dome becomes the prevalent form and 
would spread across the test surface, eventually causing an 
overheating of the test surface. 

These results were Similar to those of Patten and 
Turmeau (111, with the exception that they found vapor domes 
ceasing to appear on the surface below the 0.080 water level. 
In all levels during this study, yapor domes continued to 
feeear until forced aside by the spreading dryout without a 
dome. Although this study only went down to lnigua- de levels of 
0.040 inch, while separ and Turmeau went as low wie 0. 0207 Re 
Bpecific range was found below which just dryout without a 
dome was present. 

What causes this increase in the heat transfer coefficient 
as liquid level decreases is still not known. Pure conduction 
could cause a substantial increase in very small liqumad levels 
( = .005 inch); however, by the time this could happen entire 
dryout of the surface would most likely occur. As the liquid 


level is lowered the intense bubble-produced agitation of the 
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liguid layer decreases, and results in a lower heat transfer 
coefficient. Johnson, Jansen and Owzarski [15] suggest that 
the Marangoni effect may be the additional mechanism which, 
when combined with a bubble-produced agitation mechanism, may 
result in this increased heat transfer coefficient. Basically, 
the Marangoni effect postulates very rapid surface renewal and 
maternal circulation as the result of lateral gradzents in 
feeiace tension of the film surrounding a bubble; therefore, 
as dry spots form under a dome, Surface tension effects cause 
large velocities to exist. Additional investigation is 
necesSary utilizing visual recording techniques before any 
conclusions can be drawn as to the causes of this increase in 
heat transfer coefficient. 

When ethyl-alcohol was used as the test fluid, it was 
hard to see if dryout occurred under domes because of their 
small size and number present. In the range 0.080 - 0.120 
inches when temperatures fluctuated and the heat transfer 
coefficient increased, small domes along the edge of the test 
Surface were observed pushing fluid aside and expanding down 
to the surface itself. This same type of dome formation was 
assumed to be happening over the entire test surface, as no 
visual signs of dryout under a dome were noticed when looking 
down from the top, as in the case of distilled water. Dryout 
without a dome appeared in the 0.040 - 0.060 inch level; 
however, even with this form of dryout present nucleate 
boiling could be maintained as a result of the wettability 
Of ethyl-alcohol. Below 0.040 inch, vapor domes ceased to 
exist and dryout of the entire surface gradually occurred. 
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Be SPECIAL SURFACES TO IMPROVE THE SiEAte RA oer. 
COEFFICIENT 


ie GLhO@emcd Suurace 

A grooved surface, utilizing distilled water as the 
test fluid, was studied to determine its effect on the heat 
transfer coefficient. Results of this study are shown in 
Figure 19. An aging effect was immediately noticed between 
Pee 32 and 35. With the exception of Run 32, the effect of 
the grooved surface on the heat rae coefficient was 
Similar to that observed during tests on a mirrorlike surface 
with distilled water (Runs 10 - 22). Rue 3237S rrSst Tun 
on the new surface, has a similar profile to that suggested 
by Nishikawa, et al. [{10] for a smooth surface. 

The numerical values of the heat transfer coefficient 
are Similar to those observed earlier on a mirrorlike surface 
famede 0.2 -~ 1.0 inch liquid level. Below 0.2 inch the 
grooved surface shows an abrupt increase in the heat transfer 
coefficient, which is similar to the mirrorlike surface; how- 
ever, this increase occurs without dryout. Without the grooves, 
dryout under a dome occurred in the 0.060 - 0.200 inch range 
followed by a combination of dryout in absence of a dome and 
dryout with a dome in the 0.040 - 0.060 inch range. With 
grooves, no dryout was observed until below the 0.030 inch 
level. The observed type of dryout was dryout in absence of 
a dome. Even with this type of dryout occurring it was possible 


ee sustain mucléeate boiling down to 0.015 inchewewenedn, cue 


So 





Heat Transfer Coefficient X Lo, Btu/Hr ~ Ft or 


eZ 


jal 
_ 


[+ 
© 


=a 


eaohe: 
Run 


Rabe 


0@ Ob 


Run 


32 
eye: 
34 
Be 


} 


14.50 0g pi Eee 
10,300 Ben tee 


=/ 20 200 Suwa 


H 


=  [nGicatves, DiPyou 


e mig 
: @ © 
oe ° @ 


iG 


A 


Jali icye 


ee) 


Hrs. 


O 


Liquid Level - Inches 


O 
On , 
oO 
7p 844° 
O 
les 
Ser ea 
20 Oz 0.4 
pause ) 192 


oy 000 BU Ula 


ie 


Hrs. 


S Go = Cee a 
2 O Cun® 


Effect of a Grooved Surface on 
Heat Transfer Coefficient 


60 





eccurred on the surface, water would b@ tCranctenmecoeaien® 
the grooves to cover this dry spot. 

It appears that grooves provide a means of obtaining 
a stable high heat transfer rate at low liquid levels. No 
study was made as to the effect of liquid level on the number 
of nucleate sites or the effect of heat flux on the sites. 

No analytical solutions were carried out to relate this 
higher heat transfer coefficient to the size of the grooves, 
shape of the grooves, etc. 

2. Mesh Covered Surface 

A mesh covered surface was studied to determine its 
effect on the heat transfer coefficient at various liquid 
levels. An attempt was being made here, to Simulate a heat 
pipe to find out how the wicking material effected the heat 
transfer coefficient. The mesh, 4 layers of 40 x 50 nickel 
screen, was sintered onto the test surface as explained in 
Section III. The results of the two tests on this surface 
are shown in Figure 20. The appearance of the data indicates 
an increase on the heat transfer coefficient for depths above 
meme in going from Run 36 to Run 37. This trend conuuan 
dicts the aging data obtained previously on other test 
surfaces and is believed to be due to another mechanism. 

Run 36 was characterized by bubbles streaming forth 
between the screen layers from 7 areas along the lip circum- 
ference, followed by a stream of bubbles from the center of 
the surface once every 20 - 30 seconds. When these center 
bubbles appeared a substantial drop in the thermocouple 


temperature was noted, followed by an increase in temperature 
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until the next stream of bubbles. On the other hand Run 3/7 
was characterized by bubbles streaming forth from 9 areas on 
the lip circumference with center bubbles released every 
pe—- 3 minutes. Very minor temperaruse fluctuations were 
noted when these center bubbles were release@. The =ditterence 
in the characteristics of these two runs suggest EMat a lamge 
vapor pocket was forming within the mesh during Run 36 causing 
the surface temperature to rise substantially. ethe vapor 
pocket would then force its way through the mesh and release 
its bubbles allowing the surface to cool. On the other hand, 
the majority of the vapor in Run 37 escaped through the mesh 
wideout the sides. The abrupt increase in heat transfer at 
liguid levels below 0.2 inches during Run 36 also suggests 
some type of vapor pocket causing dryout on the surface. 

From the observations made above and the results she 
Figure 20, it appears that a substantial increase in the 
heat ‘transfer coefficient is possible utilizing a mesh ° 
covered surface, if a means for vapor escape is provided. 
As with all other surfaces tested during the study, no sub- 
stantial increase or decrease in the heat transfer coefficient 
Meeenoted in the 0.2 - 1.0 inch liquid level during any 
petveular run. No firm conclusions can be made until further 


~ 


studies on mesh covered surfaces are completed. 
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V.- CONCLUSIONS 


During nucleate pool boiling from ay horvzonea ais, 
liquid levels in the 0.2 - 1.0 inch range appear to have 
Mirtle effect on the heat eensEee COCEETC Tene: 

Below liquid levels of 0.2 inch a rapid increase in the 
heat transfer coefficient occurs with the lowering of the 
liguid level, which appears to be the result of dryout on 
mnie test surface. 

Two types of dryout occur below liquid levels of 0.2 
inch: dryout under a vapor dome and dryout in absence 

of a dome. Both types occurred with working eiaddemesd 
in this study; however, dryout under a dome was not as 
evident in ethyl-alcohol. 

Wath water as the working fluid, vapor domes are present 
in levels above 0.040 inch and domes as large as the 
boiling surface are possible depending on the combination 
of heat flux and water level. 

Aging effect present on all new or repolished surfaces 
diminishes with boiling for several hours. Higher heat 
fluxes result in faster subsiding of the aging effect. 
Reproducibility of results is possible once surface has 
elaiee. 

Grooved surfaces appear to provide a means of obtaining 


a Stable high heat transfer rate at water levels down 


eo Od Sy Suevelay 
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The scattering of data on a mesh covered surface 
indicates the presencetet vapor pockets -Wieveve.) 12a 
means iS provided for vapor escape a much improved 


performance may be possible. 


oe. 





Vl. RECOMMENDA SOR PURTHE Rs Ur 


Investigate nucleate bolting in Jueuldasle 7 -lseoeu. | 

0.2 inch utilrzing high-speed photography = eo actermmine 
the mechanism(s) which causes an increase in the heat 
transfer coefficient below 0.27 inch. 

Raise the test surface above its supports for a clearer 
visual display. 

Redesign the boiling apparatus so that the liquid level 
Can be monitored closer. This would require a larger 
imuid container to prevent the motion of the fluid during 
boiling from interfering with the liquid level 
measurement. 

Provide a fluid coolant just below he Saturation temper- 
ature to condense the working fluid's vapor. 

Mmivestigate the etrect of various gre@v7e Sizes, shapes, 
eee. On the heat transter cochiiemenitesat lower qurc 
levels. Study the mechanism(s) of heat transfer from 
grooved surfaces. 

Investigate the effect of various size mesh material, 
their shape, thickness, no. of layers, and orientation 


on the heat transfer coefficient at low liquid levels. 
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APPENDIX A 


THERMOCOUPLE CALIBRATION PROCEDURE 


Six metallic-sheathed copper-constantan thermocouples 
were utilized to measure the temperature gradient in the 
nickel cylinder and ultimately, through the use of the least 
Squares method, to determine the temperature at the boiling 
surface. The properties of these particular metals proved 
them to be most suitable for the temperature range studied. 
A total of 8 thermocouples were calibrated to take care of 
any potential thermocouple problems. 

To insure accuracy in meaSurement, the thermocouples 
were calibrated at fixed reference points spanning the 
experimental temperature range. These fixed reference 


points as recommended by Roeser and Lanburger [16] included: 


Mm@petriple point Of ice -—= O20leCe (3270275) 

2) Boiling point of water -- Dependent upon atmospheric 
pressure 

3) Freezing point of -—pure: tins=—» 23198826 2(449..3828) 

4) Freezing point of pure lead -- 327.4°C (621.32°F) 


All calibration voltages were measured with the reference 
junction in a distilled water-ice bath and read with a Hewlett- 
Packard Data Acquisition System Model 2010C containing a 
guarded data amplifier and an integrating digital voltmeter. 
This instrument gave a printed digital output with an accuracy 


Siete 0.5 MICLOVOLES. 
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To measure the triple point of ice an ice bath, con 
Sisting of finely crushed ice made from distilled water, 
was contained in a 16-ounce thermos bottle and saturated with 
distilled water. Three determinations of the voltages for 
each thermocouple were made and the results averaged and 
iliecad in Table 2. 

To measure the boiling point of water, a stoppered 600 
milliliter glass beaker, vented to the atmosphere and filled 
with 300 ml distilled water, was heated with a standard hot 
plate until rapid boiling occurred. The water was allowed 
to boil for 15 minutes to remove entrapped air. The 
thermocouples were then inserted individually through the 
stopper to a point three-fourths of an inch above the boil- 
ing water surface, and at this point a valid reading of the 
vapor temperature voltage was measured. Three readings of 
voltages at this particular pressure were made and the results 
averaged and listed in Table 2. 

To calculate the equilibrium temperature of boiling 
water the Fees phere pressure was recorded from an aneroid 


barometer and the following equation used [16]: 


3 
To,- = 100 + 28.012 (P - 1) - 11.64 (P - 1)* + 7.1 (B - 1) 
12 134 12. 
O O Oo 
where 
eat = equilibrium (Saturation) temperature in degrees C 
P = local atmospheric pressure in mm Hg 
ES = standard atmospheric pressure in mm Hg 
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To measure the freezing point of tin, a National Bureau 
of Standards tin sample {melting peomnt — 231 oe ee ac 
heated in a stoppered test tube in a Hoskins Electric 
Furnace. The thermocouple was inserted through the stopper 
in a glass protective tube into the molten tin to a depth 
of approximately 1 1/2 inches. Powdered graphite was sprin- 
kled on the molten tin to retard oxidation. Once the tin was 
thoroughly molten, the furnace was secured, and the test tube 
removed and placed in an asbestos lined container. The 
temperature of the molten tin in terms of thermal emf of 
the thermocouple was then monitored with the Data Acquisition 
System until the freezing point was reached and passed. The 
freezing point of tin was characterized by an almost constant 
voltage (read every second) for a period of approximately 2 
minutes. As noted by Von Vlack [17] the temperature of tin 
dropped below the freezing point and returned to a constant 
temperature which was considered the freezing point. Once 
the freezing point was determined, the test tube was placed 
back in the furnace and heated for a second time and the same 
process was carried out for another temperature reading on the 
same thermocouple. Each thermocouple was read twice and the 


results averaged and listed in Table 2. 


In a similar manner the freezing point of lead was 
measured using a National Bureau of Standards lead sample 
(melting point 327.4°C). However, with lead no constant 
voltage was observed for longer than 10 seconds, therefore 


a millivolt vs. time curve was plotted and the inflection 
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point found and an uncertainty assumed around this point. 
A second set of data verified this point and the point with 
ims uncertainty Jlistedmin taplesz- 

Using the voltages obtained for the four fixed reference 
points to solve for the four unknown constants (Appendix B), 
the following equation was used to determine a temperature- 
voltage relation for each thermocouple: 

T = A+ BM + CM* + DM? 

where 


db thermocouple temperature in degrees F 


M 


thermocouple voltage in millivolts - mv 
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APPENDIX B 


SAMPLE CALCULATIONS 


1. Thermocouple Calibration (thermocouple #20) 

Utilizing a Least Squares Polynomial Fitting 
Smroutine (LSOPL2) for an IBM 360 computer, coefficients 
for the power series: 


24 px? 


T(X) = A + BX + CX 
were determined from obServed dependent and independent 
waiereroles found during the calibration of the thermocouples. 
For thermocouple 20 the observed dependent and independent 
variables were respectively: 

a) The fixed calibration temperatures in degrees F; 

S(O: 
22 o4 
449.38 
G2. 32 

b) Millivolts from calibration process corresponding 

wOuthesnuxed calibration. eemmeracunes: 
eel 

4.314 
One 9 
16.450 

With these variables the resultant power series for tempera- 
eewe (T) as a function of millivolts (xX) for thermocouple 
20 was: 


T(X) = 30.12597 + 45.74406X - 0.88978x? + 0.01786x° 


Similar solutions were obtained for the other thermocouples. 
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2. Data Reduction Program 


With the aid of a data reduction program written for the 


IBM 360 the following values were determined: 


sa - equilibrium (saturation) temperature, °F 

ae - surface temperature, 3B 

TT. ~ Tsatm surface-Saturation temperature difference , °P 
k = Conductivity, Btu/Hr —9ee °F 

(Q/A) - heat flux from power input, Btu/Hr - Ft? 
(Q/A) in - heat flux from temperature profile, 


Btu/Hr - Ft? 
=) , was determined from the equation: 


Sa 


3 
T=) C0) 205012 m= oA (er = ae) aes lie eel 


Po De Po 


ex. Run 16 where P = 761.90 mm Hg 


To., = 100 + 28.012(761.90 - 1) -11.64 (761.90 - 1) 
760.00 7E0n 00 
ASG Ome 
760.00 
fee = 100-07°C (212.12°r) 


b) T,, was determined by utilization of the method of 
least squares for solving the equation: 
i aX + B (see Appendix D) 
CxmwhUn 6) — DeieaPointe & 
With the necessary information from smooth data, 


the following sums were formed: 


Ue 





1 V8 
262.16 0.220) 
251.46 0.715 
238.99 lees 
Se = 752.61 SS ro 
TAG x54 
57.68 .048 
179.79 Poul 
200.37 > lea a6 
pee —=527.04 SX. 7 = 922005 
i 4 1 


SuDSt LEViEInGmaMEOmENeweaguatteon fom and £: 


ie 3 =X.T, -~ ( 2X.) ( ele) 
3 apne Meese yo) a 
BE le 
Pen 3 527.84) = (2 CSO S256 1 ea 5 
a (2 4025) -— "9 (42622) 1.483 
Boge (Soar 0 ok a = (29 nt Sx.) 
7 nt a iii i 
Z 2 
Be= (752.601 1420035) 15277 e4 2 oO SoG = ome 


B (2 2065) aed 622 1.483 
Substituting into the original equation: 
LS aS ae HST), oi0 
and solving for the temperature at the surface (X = 1.710) 


Tt = 22). 62go 
W 
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Sh determined from: ex. Run 16 - Data Point 3 


eee Psat 
227.627—- 212 ec Ona 

qd) Thermal Conductivity - An equation for the thermal 

conductivity as a function of the temperature was 

obtained from the following data for Nickel 200 

meormerly Nickel “A™) (He) utilizineg@agarinme tie 


Least Squares Polynomial Fitting Subroutine for a power 


series: 
= 2 
K(T..,) = teas Rea 
AGB) = 16 Nenalstly eo syersie icles sls eielaljolspaeheh ate: Al 


Thermal Conductivity of Nickel "A" 


TOF k (Btu ) 
Hr - Ft - °F 
BO Oe A 
UNG B7 5 
392 33.4 
592 29.0 


The coefficients obtained resulted in the following 


equation: 


2 


ns de, SO SSS 4 a AO ACIS AM oo -O000000077T_., 


Where se is the average temperature over the interval 
(T, + T, + T ) 


3av 
5 


ex. Run 16 - Data Point 3 (ie =P Once) 


k(T ) 


~ 42.8692 =" ,0241952( 250337) 0000000077 (250.87) 


KAT oy) 36.30 Btu/Hr - Ft - °F 


ies 





See 


e) Calculations of Heat Flux from Power sino: 


(2 = Voltage xX Ciirrent < 3.43 eee. 4 
A), Hr - WATT prt? 


ex. Run 16 - Data Point 3 - 27.5 Volts, 3.20 Amps 


(Q) = 27.5 X 3.20 X 3.413 = 15,807.57 Btu/Hr - Ft? 
A 


Oe, O19 


ft) Calculation of Heat BPlux from Temperature Profile; 


ex. Run 16 —- Data Point 3 
(2) = 36.80(262.16 - 238.99) °F x 12in = 10,283.29 Btu/Hr-Ft* 
A (Hr-=Pte °F) xX .995 an xX Ft 


Heat Transfer Coefficient was Determined from the Equation: 


b= (0/A) nose 


a - weer 


where (0/A)m_.,, = average heat flux determined from the 
BeMDekaLuLe proLile GquEing salnececmn— 
ikeleell. ie bly 


nh 
(2 = (Q/A) ,. eye 
- 1 


n = No. of data points in a given run. 


es. Run 14 —- Data Point 1 


0 = 14,200 
Cz) 
T-av 
AT = TT - t= = 22.91 
aS IL Ge = 620 Btu/Hr - aes > |p 
js S) AL 
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APPENDIXME 


UNCERTAINTY ANALYSIS 


In order to check the reliability of the data collected 
and to determine the magnitude of the most probable error an 
uncertainty analysis was performed utilizing the method of 
Kline and McClintock [19]. 

A. Uncertainty of Measuring Distances: 

ir. | Dhesnocousre poOsctelens = =o. 00045 t 


AX = Ay = X3 


2 Zhe: = 
AX oe + ne = 2 X .00000016 = .00057Ft 


2. Liquid Level: 


W 


Low Heat Flux-- Q/A 215,000 Btu/Hr-Ft? 
Wrp- t OG) Ole t 
2 
Medium-Low Heat Flux-- 15,000 Beu/He-Fe- ZL. Oye 25 000 Btu, Bete 
Wo mee) OUP Zea 


2 


Medium Heat Flux-- 25,000 Btu/Hr-Ft~ v O/A v_ 60,000 Btu/Hr-Ft? 


Wor TOOSOS Et 
B. Uncertainty in determining heat flux from temperature 


profile: 


Q=k AT, 
AX 


where k = conductivity, Btu/Hr-Ft °F, 
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Uncertainty in temperature from calibration curves 
ise+ 0.2°C fromTablle 9 ineRetexeneami 16). 


Wr = “fa. JO 2e 


UNCErtainey snes. =the average temperature 


3av 


of the 4 thermocouples located in position X3- 





Tag fa bee 
4 
W = Ww = W =< Vy = + .36 °F 
le Ty dts Ta) 
OT3ay =7 3av OT 3 ay ol 35y 2 
W — 
Tay 





T3av a 


“5 al, ST ae) ee emer 
av a 


Uncertainty in (TS = Dag 








ATy ~ oa i‘ Toes 
oAT By Ass 

ee Xx = —] 

Ons) oT 304 
2 2 2 2 

ae aa (W, ) + (-w a (.36) + (-.18) 

Xx 1 T3av 

= + .40 °F 


78 





Uncertainty in average temperature of the Nickel 
cylinder for determining conductawity: 


At Tiamat err A 





av i Z 3av 

ON 2 olay = olay ee 

oly Oe OT 35v 3 
W = WwW = 2w 

Th, T'5 lv 3av 


= / 2(.36)* + (18)? 


av L T3av 3 


= 
= 
I 
WOIND 
3 
Ke 
a — nN) 
-\- 
Wie 
= 
NO 


Uncertainty in determining condwerwmavity—— Kk: 


k(T,.) = 42.8652 - .0241952 T_ + .0000000077 (T__)? 
av av 
ok = —~.,0241952 + 2(.,000000007taaT > 
Ai av 
BY) 
ex. Run 16 Data Point #3; ‘ee = 250.87 
ok = —-.0241952 + 2(.0000000077} (250.87) 
olay 
= 
We — ( ae Wr \4 Se ae vin 
3 av ay OT aw av 
=m (O24). 18)y = + .004 Btu/Hr-Ft °F 


The above uncertainty in conductivity was determined from 


the equation for bi OUBe above. The data from which this 


equation was formed, using the least squares polynomial curve 


fitting technique, was believed to have an uncertainty of + 1% 


which is considerably larger than the uncertainty given above; 


therefore, the actual Uuncertainey in sae was assumed to be: 
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Uncertainty in pressure readings from aneroid 
barometer; 


Wo = + -5 mm Hg 


Uncertainty in determining saturation temperature: 
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D. Uncertainty in heat transfer coefficient: 
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Uncertainty at a higher heat flux: 


ex. Rune! 6 Data Pernt 5 
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AT = 23.93 °F 
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av 
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APPENDIX D 


LEAST SQUARE METHOD FOR OBTAINING SURFACE TEMPERATURE 


Figure 21 is a plot of the temperature distribution 

within the nickel test cylinder at various power inputs. As 
shown the temperature profile appears to be linear, and the 
surface temperature (temperature at point 0.0) could easily 

be determined by extending a line through the plotted points 
to the surface; however, this would have been a time-consuming 
process and another method was sieghwalnl The linearity of 
the temperature profile provided an opportunity to obtain a 
functional relationship for the surface temperature by using 
the method of least squares as described by Holman [20]. An 


equation of the form: 
T = AX + B 


was required to obtain ees where the constants A and B were 


determined from the equations: 





ae Sea OF 5 ee Lae 3 
3xx* - ( =x.)? 
al at 
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( a | = X;) - ( wo = Xj) 


2 2 
3=X, - ( =X,) 
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With the axis of the coordinate System at the bottom of the 


nickel cylinder, the following values were determined by the 


measurement of thermocouple position, and surface height. 


Lest, Cylinder s+ Test Cyijsmder 42 
X, == 15 ara) aig). Xy = 0. A@F in. 
Go Og Tbe ea iele 45 = Osean: 
oS tine xa = Lf 20e 4 nie 
Ae leet LO) inne X, = 1.690 in. 


With the constants determined Gite 


a 


profile and the distance to the surface 


given temperature 


known, the surface 


temperature, Wee was easily determined as shown in Appendix 


Be 
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